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ABSTRACT 

The size and cost of microcomputers continue to decrease 
while their memory capacity and execution speed increase. 
These advances should result in small, inexpensive machines 
attaining the same computing power as current mainframe 
models. The interim need is to adapt general finite element 
codes to present day, less capable microcomputers. This 
thesis explores the program structure, memory management, 

I/O procedures and equation solving methods necessary to 
accomplish that task. The equation solving capacity and 
speed of the Apple-II Plus Personal Computer Systems and the 
Hewlett-Packard System 45(A) Desktop Computer are compared. 

A finite element program for the static analysis of space 
trusses is presented, as adapted to and tested on the Apple-II 
Plus. The program output may be printed in either English or 
French. 
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I . INTRODUCTION 

A. PURPOSE AND SCOPE OF THE INVESTIGATION 

In the preface of the third edition of his classic refer- 
ence on the finite element method, Zienkiewicz [1] made an 
assessment of the "state-of-the-art” of this most useful analy- 
tical technique. In 1977, Zienkiewicz described a discipline 
which had ammassed a collection of nearly 8000 published ref- 
erence works in an ever widening variety of applications, 
and which had evolved far beyond its roots as an engineering 
tool for linear analysis of solid mechanics structures. 

Present applications commonly include heat transfer, fluid 
mechanics and non-linear solid mechanics, as well as the com- 
plete static and dynamic analysis of engineering structures. 

A more geneal statement may be made [1] in describing the 
finite element method as " . . . a general discretization pro- 
cedure of continuum problems, posed by mathematically defined 
statements." The method now becomes the means of approximat- 
ing a continuous problem with a discrete model composed of 
finite elements, as opposed to "infinitesimal" elements in 
the sense of the calculus. In some cases, in the limit, as 
the number of elements increases to infinity and the physical 
dimensions these elements represent decrease, the solution 
of the model becomes exact. It is this generality of method 



11 



and application which has been one of the two most important 
factors influencing the widespread use of the method in 
scientific research and industrial applications. 

The second factor influencing the spread of the finite 
element method has been the availability of high speed, high 
capacity, general purpose digital computers. It is well 
known that the method remained in its infancy until general 
purpose digital computers became available in reasonable 
numbers. Today technological advances in microelectronics 
and solid state devices have drastically reduced the real 
cost of computer time. Advances in the art of building 
static, random access memory now allow the use of Direct 
Memory Access techniques in small (16+ bit word) desktop 
machines. When these machines are interfaced with magnetic 
hard disk drives, a small user easily has the potential of 
16 Megabytes (8 binary bits equals one byte) of high speed 
storage, at a small fraction of the cost of a large main- 
frame computer [2]. These developments are truely revolu- 
tionary. Very "friendly" and general structural finite 
element codes, some with appealing interactive graphics 
capabilities, may soon become available to the smallest of 
o r gan iza t ions . 

In addition to the new machines now on the horizon, there 
are a very large number of smaller, less capable computers 
presently in service. These smaller machines, often des- 
cribed as "microcomputers," vary in word length and storage 
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capacity, generally having from 16 to 256 thousand bytes of 
random access memory. They support a wide variety of periph- 
eral equipments, the most ind isp ens ib le of which are high 
speed, magnetic disk drives . ^ These smaller machines are pro- 
portionately less expensive than larger computers, which 
accounts for the very great demand for these systems in industry, 
government and the private sector. Many of the most popular 
models began as "personal’ 1 or leisure computers for hobbiests, 
but rapidly moved into small business, medical and financial 
management applications as the utility of the machines was 
realized and software became available. 

If one looks, however, for finite element software or 
attempts to find evidence of this widespread analytical tech- 
nique being implemented in any substantial manner on micro- 
computers, the results are disappointing. It is the thesis 
of this investigation that: 

(1) Substantial problems in finite element analysis can 
be solved on microcomputers. 

(2) Microcomputers can solve substantial finite element 
problems with acceptable accuracy and to reasonable precision. 



Disk storage extends the usefulness of microcomputers by 
avoiding or softening the impact of some of the limits 
imposed by having too little main "core" memory. 
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(3) Although the central processor clock speed of most 
popular microcomputers is much less than that of a mainframe 
computer, and in spite of the interpretive nature of most 
microcomputer languages, substantial finite element problems 
can be solved on these machines in an acceptable amount of 
time. 

(4) General guidelines can be established for using inno- 
vative program structuring, memory management, software overlay 
schemes, data base management and high speed, disk storage to 
allow the implementation of substantial finite element codes 
on small popular microcomputers. 

3. CHOICE OF MACHINES AND PERIPHERAL DEVICES 

The choice of machines was influenced most heavily by the 
types of systems available at the Naval Postgraduate School. 

A Hewlett-Packard System 45 Desktop Computer was available 
within the Mechanical Engineering Department and was a logical 
first choice. Previous work had already been done on imple- 
menting a finite element code on this system [3] and the 
potential existed for increasing the capabilities of the 
program by adding additional finite elements, changing from 
an "in-core" to an "out-of-core" equation solver, and using 
magnetic floppy disk instead of magnetic cartridge tape mass 
storage. These systems are not uncommon within the Department 
of Defense and the Navy and the potential benefit from an 
improved program was judged to be worth the additional effort 
which would be required to upgrade the code. 
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1 he System 45 is, however, nearly outside the boundaries 
of the group of machines previously described as "popular 
microcomputers," The base price of the computer alone, ex- 
clusive of the disk drives, was over $29,000.00 in 1977. A 
new System 45 (HP 9845C), upgraded and having slightly more 
computing capability, would have cost about $39,500.00 if 
purchased in 1980 [4] and still be without disk mass storage. 

It was decided, therefore, to also choose a system which 

would represent those on the lower end of the cost spectrum. 

2 

Several Apple-II Plus' Personal Computer Systems were avail- 
able, having a total price, including multiple disk drives, 
of approximately $6,000.00. In addition to having a low cost, 
the Apple-II Plus systems also had the capability of the 
FORTRAN Language Option, not available on the HP 9345A. 

For the reasons enumerated above, it was decided to in- 
vestigate the equation solving capability of the Apple-II 
Plus Personal Computer and the Hewlett-Packard System 45 
Desktop Computer (HP 9845A) and to attempt, should adequate 
time be available, the installation of representative finite 
element codes on both systems. 

1 . The Hewlett Packard 9845 Desktop Computer 

The System 45 was, in 1980, the t op- o f - t he- 1 in e model 
of Hewlett-Packard desktop computers. As with all of this 



“Apple-II Plus is a registered trademark of the Apple 
Computer Company. 
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firm s computers, technical details about the hardware 



are 



closely held by the company and few are available in users 
manuals supplied with the machine or in the Hewlett-Packard 
Jou rna 1 . References [5] and [6] and the installation docu- 
ments for the system provided the information necessary to 
conduct the investigation. 

Only a very brief summary of the most important details 
appears in the following sections, and readers intending to use 
programs presented in this thesis must become familiar with 
the machine in general and their own systems in particular. 

a. System Configuration 

The system configuration during the investigation 
is shown in the table below. Only the minimum devices and 
options which are required to run the program presented in 
this thesis are shown. Note that the possession of a more 
capable system than the one shown (e.g. one hard disk instead 
of dual flexible disk drives) will not preclude use of the 
proposed software. The reader need only change the "mass 
storage unit specifier” [6] and "select" [5] codes to match 
his/her particular configuration. 

Although, under the "unified mass storage con- 
cept” [5], it is possible for a user to substitute cartridge 
tape for disk mass storage in nearly any program, this prac- 
tice is unsuitable for solution of large systems of linear 
equations out-of-core. Previous experience with less demand- 
ing mass storage tasks [4] showed degradation of system 
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performance, on large problems, where tapes were required to 
rewind too many times. Stretching of the tape initially 
introduced parity errors and later caused breakage of the 
tape. 



Table 1: Configuration of the Hewlett-Packard 

Systems 45 Installation 



Required Devices: 
Part No. 



Des cr ipt ion 



Remarks 



HP 9845A 
HP 9 80 32 A 
HP 9885M 
HP 9 8 8 5 S 



Computer with CRT 

16 Bit Interface (9885A Disk) Option 85 
Flexible Disk Drive (Master) Drive 0 

Flexible Disk Drive (Slave) Drive 1 



Required Options: 
Numbe r 



Descript ion 



203 


64K RAM 


310 


Mass Storage ROM 


320 


1/0 ROM Left 


330 


1/0 ROM Right 


370 


Graphics ROM 


500 


Internal Thermal 



Remarks 



Select Cede 13 
Select Code 0 



Select Code Summary: 

Code Device 

0 Printer 

8 Flexible Disk Drive (Controller) 

13 Graphics ROM 

16 CRT Display 
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Machine Precision 



b . 

Tne HP 9 84 5 A stores floating point numbers between 
9.99999999999 X 10 and 1 X 10 in absolute value. An ex- 
tended calculation range of 9.99999999999 X 10 511 to 1 X 10~ 51 
is used for the results of intermediate calculations [5]. 
Floating point numbers are described as full precision or 
short precision. Full precision floating point variables 
each occupy 8 bytes of storage and are used exclusively in 
all HP Enhanced BASIC Programs presented in this thesis. 

Chapter 3 of Reference [7] discussed the estima- 
tion of machine roundoff error (as part of a discussion on 
estimating the relative error in the solution of systems of 
linear equations from the matrix condition number). The 
worst case roundoff error, in base 10 for a 12 digit machine 
is : 

10~ 12+1 = 0.00000000001 

Section 2.2 of Reference [7] also defines a 
quantity called the machine epsilon as the smallest float- 
ing point number distinguishable from zero. A program, 
modeled after that presented in the reference, found the 
machine epsilon of the HP9845 to be 7 . 2 7 5 9 5 7 6 14 1 5E - 12 . 
c. The Enhanced BASIC Language 

The Enhanced BASIC Language from Hewlett-Packard 
has all the features of the most powerful, high level versions 
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of BASIC (Beginners All-purpose Symbolic Instruction Code) 
and certain extensions designed to provide optimum use of 
the System 45. References [5] and [6] contain a complete 
description of the language and references for further study 

As with most versions of microcomputer BASIC, 
the language is completely interpretive in nature. That is, 
no compilation or conversion of source code to machine lan- 
guage is accomplished before executing the program. Each 
line of source code is interpreted, line-by-line, as execu- 
tion progresses. Although this makes it possible for the 
program to have such incredible flexibility as to be able, 
in some cases, to edit itself during execution, it also con- 
siderably slows the speed of execution. On the positive 
side, it is true that there are many instructions in Hewlett 
Packard Enhanced BASIC which can reduce the amount of source 
code required over, say, the same program in FORTRAN (e.g., 
the MAT series of instructions for matrix manipulation [5]). 

The reader who is familiar with FORTRAN will 
have little trouble following the programs presented in this 
thesis. The table below lists some of the more unusual fea- 
tures of the language and should be enough to allow complete 
understanding of program flow. For further information, the 
reader must consult the references given above. 
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Table 2: Selected Members of the Hewlett-Packard 

Enhanced BASIC Instruction Set 



Def init ions : 

"f ilenumber' 1 is a numeral or numeric expression 
having an integer value between 1 and 10. 

"file specifier" is a character string enclosed in 
quotation marks or a string variable which is a valid 
and complete file name (including mass storage unit 
specif ier) . 

S ta t emen t s : 

ASSIGN // filenumber TO file specifier 

The ASSIGN statement places the filename in an 
internal table and allows the file to be referenced by 
a single integer number in PRINT// or READ// statements 
[5] . 

BUFFER// f il enumb e r 

Sets up an additional 256 byte I/O buffer for all 
READ// and PRINT// statements to the specified filenumber, 
allowing the most efficient coupling between a file and 
its associated mass storage device [6]. 

OPTION BASE 1 

Specifies that the lower bound of all array variable 
subscripts shall be 1 (e.g. Vector (1) is allowed but 

Vector (0) is not) [5]. 

Yyy = FNXxxxxx (aa, B,Ccc, . ..) 

Call to a function subprogram previously defined by 
a DEF FNXxxxxx statement [5]. In FORTRAN this would be 
YYY = XXXXXX (AA, 3 , CCC, . . . ) where XXXXXX is a FUNCTION 
subprogram or an intrinsic function set up with a 
statement function definition. 

ZzzOOO : 

An example of a label to which program flow may be 
transferred and which keeps its relative position within 
the file regardless of changes in line numbering [5]. 

Any proper statement may follow the colon. 

BEEP 

Causes a tone to be sounded at the keyboard [5]. 
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DISP "Character string" 

Causes the character string enclosed in quotation 
marks to be printed in the display area of the CRT [5]. 

OVERLAP 

Allows computation and I/O to occur simultaneously. 
Most beneficial when a program has nearly equal amounts 
of computation and I/O. This mode is disabled by the 
SERIAL statement. See [5] and its references for a 
discussion of the OVERLAP statements impact on error 
traps . 



When attempting to reduce the run time of a pro- 
gram, it is important to recall the method used by the HP9845 
when a subroutine or function call, or transfer to a state- 
ment label or line number, is encountered. The BASIC inter- 
preter jumps to the beginning of the file whenever a search 
for one of these items is instituted. Run speed can be im- 
proved significantly when the most frequently used functions 
and subroutines are placed as near to the beginning of a long 
program as possible. Note also that comment lines in BASIC 
(statements preceeded by REM or an exclamation point) are now 
a liability because they increase the length of a program. 

The amount of documentary comments should be limited to the 
minimum necessary to understand the code. 

It is possible to overlay subroutines, by using 
the LINK statement, at any line number in the main program 
(or in place of the main program if desired) [5]. Overlaying 
reduces the amount of program in core, allowing the storage 
of more variables or the execution of programs too large to 
fit in memory :_n a single piece. 
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d. Disk Mass Storage Considerations 

The characteristics of the disk mass storage 
hardware used with any microcomputer have a significant 
impact on the performance of the system in finite element 
calculations. The size and number of usable physical records 
on the disk, size of the file directory available and types 
of file structure and access allowed are all important con- 
siderations in attaining optimum use of out of core storage. 

In the HP 9845 there are two types of records with 

which the programmer must be concerned [6]. A phy s i ca 1 record 

* 

is the unit of storage dealt with by the mass storage device 
itself (and is usually the buffer size associated with the 
device). Def ined records are the smallest addressable units 
of storage, in even numbers of bytes, which can be individually 
accessed by the user. Whenever possible the defined record 
length should match the physical record length of the device 
in use. In the H? 9345 this matching of record lengths can 
provide an improvement in I/O performance of two- or three- 
t o-one . 

Reference [6] describes the use of binary DATA 
files and Direct Memory Access (data transfer without the 
use of buffers) for the rapid transfer of entire arrays 
between disk mass storage and core memory. This method can 
result in a considerable I/O improvement over buffering tech- 
niques, depending on the storage device and the amount of 
data to be accomplished. 
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The characteristics of the HP 9885 Flexible Disk 



Drive are listed in the table below. 



Table 3: Selected Characteristics of the 

Hewlett-Packard 9885 Flexible Disk Drive 



Storage Capacity: 

Bytes (total)- ---------------- 499,200 

Physical records (maximum) ---------- 1,950 

Physical records (useable) ---------- 1,901 

Files (Director size limit)- --------- 352 

Bytes per physical record ---------- 256 

Accessing: 

Maximum transfer rate (bytes per second) - - - 46,000 



2 . The Apple-II Plus Personal Computer 

The Apple-II Plus Personal Computer is manufactured 

by Apple Computer Incorporated of Cupertino, California and 

is designed around the MOS Technology 6502 microprocessor. 

1 6 

The maximum addressing range is 2 (64 k) locations on a 

16 bit parallel bus. The data bus is S bits, parallel and 
bi-directional and the CPU clock speed is 1.203 MHz. More 
technical information is available in Reference [8]. 
a. System Configuration 

The baseline Apple-II Plus central processor 
does not vary appreciably from installation to installation, 
although enterprising companies have recently marketed acces- 
sories which can actually replace the 6502 microprocessor 
with the more popular 280 CPU (and/or increase the clock 
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speed). The possible systems, however, and the languages 
available on those systems are too numerous to mention. One 
special purpose and six general purpose "slots” are available 
at the rear of the computer for connection of peripheral 
devices. These slots are supplied power and have access to 
all data bus, address bus, control and interrupt lines. 

The table below describes the configuration of 
the system used in this investigation. 



Table 4; Configuration of the Apple-II Plus 
Personal Computer Installation 



Ha r dwa re : 

No. used Device 

1 Apple-II Plus 

1 Apple Language System 

2 Apple Disk-II Interface Card 

4 Apple Disk-II Floppy Disk Drives 

1 Apple Parallel Printer Interface Card 

1 Microline-80 Dot Matrix Graphics Printer 
1 Sony Trinitron Color Video Monitor 



Remarks 
48K RAM 
16K RAM 



Software: 

Apple FORTRAN 
Apple Pascal 



Slot Connections: 

Slot No. Device 

0 Language System Card 

1 Printer Interface Card 

5 Disk-II Interface Card 

6 Disk-II Interface Card 



Rema r ks 



Volumes 11 and 12 
Volumes 4 and 5 



b. Machine Precision 
The Apple-II Plus 

j_ 3 s 

between 1.7 X 10 and 5.3 X 



stores floating point numbers 
— 3 8 

10 in absolute value and 



integer numbers between +32,768 and -32,768 [8]. Floating 
point variables each occupy 4 bytes (two machine words) of 
storage and integer variables occupy 2 bytes (one word) . 

The bit allocation for a floating point number is: 

31 30 ... 23 22 ... 16 and 15 ... 0 

Item: sign exponent mantissa 

Floating point variables having a maximum of 6 decimal digits 
may be stored. 

Using the methods of chapter 3 of reference [ 7 ] , the 
estimate of worst case roundoff error is: 

10~ 6+1 = 0.00001 

The machine epsilon (smallest floating point num- 
ber distinguishable from zero) was investigated using the 
program of section 2.2 of reference [7] and was found to be 
. 119209E-06 . 

c. The Apple FORTRAN Language and the Pseudo Machine 

3 

References [9] and [10] discuss the Apple FORTRAN 
language and explain its implementation on the machine and 
relation to the Apple Pascal Operating System. Although 
it is not necessary to know any Pascal to operate the system 
and use Apple FORTRAN, from a user point of view it is help- 
ful to remember that this version of FORTRAN was implemented 

3 

Apple FORTRAN is a registered trademark of the Apple 
Computer Company. 



25 



